Entropy and light cluster production in heavy-ion collisions at intermediate energies 

Yogesh K. Vermani and Rajeev K. Purij 

Department of Physics, Panjab University, Chandigarh -160 014, India 

(Dated: July 20, 2010) 

The entropy production in medium energy heavy-ion collisions is analyzed in terms of ratio of 
deuteronlike to protonlike clusters (diike/puke) using quantum molecular dynamics (QMD) model. 
The yield ratios of deuteronlike-to-protonlike clusters calculated as a function of participant proton 
multiplicity closely agree with experimental trends. Our model predictions indicate that full ther- 
modynamical equilibrium may not be there even for the central geometry. The apparent entropy 
extracted from the yield ratios of deuteronlike-to-protonlike clusters, however, reflects the universal- 
ity characteristics i. e. it is governed by the volume of reaction independent of the target-projectile 
f^ ' combination. Our calculations for apparent entropy produced in central collisions of Ca-|-Ca and 

vN ' Nb-|-Nb at different bombarding energies are in good agreement with 47r Plastic Ball data. 



I. INTRODUCTION 

Various experimental and theoretical studies in the recent past have indicated a clear demarkation of colliding 

matter into participant and spectator matter especially at high incident energies [l|, l2|. This is characterized by 

^ ' the formation of hot & dense fireball. In addition, the size of the participant volume or 'source size' is also linked 

'Zlf '■ n n 

^ , with the emission of composite particles and ultimately vifith the mechanism of the production of entropy 3-5]. The 

highly dense fireball formed for a very short duration (< lO^^^s) can be of importance for probing the properties 

►-j^ of condensed nuclear matter. The formation of the fireball is reported to be affected by many factors such as beam 

l/^ , energy, overlapping volume, as well as density reached in a reaction. Earlier experimental studies using the Plastic 

^j^ , Ball/Wall detector [y, |7| showed that the deuteron-to-proton yield ratio varies with the impact parameter of the 

reaction indicating a strong dependence on the participant volume. Such behavior is believed to be resulted from the 



C ■ coalescence production mechanism [6|, |7| . Experimentally, the study of d/p ratios obtained in asymmetric collisions of 

^^ P+Kr, 0+Kr, and Ne+Ar also indicates the importance of coalescence mechanism [8[. In this study, it was observed 

^-H . that d/p ratio increases with the beam energy. The baryonic entropy, however couldn't be directly measured from 

I> . these ratios, since total yield is required from the source. The measurement of the yield ratios and baryonic entropy 

K> , produced are, therefore, very promising candidates for estimating the fireball produced in the hot & dense nuclear 

H , matter. 

As a matter of fact, the pion production also contributes significantly |9l4l2| towards the entropy generation at SPS 



[ISlllJI and higher energies. At SPS energies (^^160 AGeV fixed target), the pion number increases with beam energy 
to about ten times the number of the nucleons [IJ]. There, an increase in the entropy production was observed with 
beam energy as one moves from AGS energies towards SPS and higher energies [10|. This enhancement may also be 
conjectured as manifestation of the change in the collision dynamics at such high energy 15|, |l6|. For the incident 
energy range considered (400-1050 AMeV ) in the present work, the inclusion of the pion production is not going to 



m the pre 



affect the entropy production appreciably [Ul, ll2|, [17[ . For instance, at SIS energies (upto 2 AGeV), the total number 
of pions is only 10 % of the nucleons [18]. 

In the present work, we aim to estimate the baryonic entropy Sn generated in the fireball via the yields of free par- 
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tides and light clusters {i.e. protons, neutrons, deuteron, tritium, helium-S, and a-particles). We made a systematic 
study of symmetric heavy-ion (HI) collisions in the incident energy range 400-1050 AMeV within dynamical model 
namely, quantum molecular dynam,ics (QMD) model [19| . We shall also compare our results with experimental data 
taken with Plastic Ball/Wall detector. The dynamical approach such as QMD model provides an useful platform to 
estimate the yield of composite particles in terms of n-n correlations from the start when the colliding nuclei are well 
separated to the final state where matter is cold and fragmented. As reported in Ref. [20] , the production of light 
charged particles carries vital information about the stopping and equilibration of nuclear matter during a collision. 

II. THE QUANTUM MOLECULAR DYNAMICS MODEL 

The quantum molecular dynamics (QMD) model is a n-body transport theory that describes HI reactions in 
intermediate energy regime (20 MeV/nucleon < Eiab < 2 GeV/nucleon) on event by event basis. The two essential 
ingredients of the model are nucleon-nucleon (n-n) potential and stochastic scattering. Here each nucleon follows the 
trajectory according to the classical equations of motion [19[ : 

a(H) dm .,. 

The expectation value of the total Hamiltonian ("H) consists of kinetic and potential energies terms as: 

in) = (r) + (v) 



= ^^-f^(^/'^«™^ + v;^«'= + yf°«' + \//'^^^). (2) 

2 — 1 2—1 

The 1^ yrme^ V^^'^, vf^°^^ and V^^'^^^ in Eq.(I2|) are, respectively, the local Skyrme, long range Yukawa, an effective 
charge Coulomb and momentum dependent parts of the interaction |l9|. The potential part without momentum 
dependent interactions (MDI) is called static equation of state (EoS). For the present study, we employ a soft equation 



of state represented by incompressibility k = 200 MeV and energy dependent n-n cross section [19j. The equation of 
state with MDI is labeled as 'SM' whereas without MDI, it is labeled as 'S'. Since many studies reveal the nuclear 
matter to be softer 
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22| . we also prefer the soft EoS. 
The entropy information may be obtained from a classical charge symmetric oas of nucleons and deuterons in 
thermal and chemical equilibrium using the relation suggested by Kapusta et al. [3|, l23| : 



■>N 



3.945 - InRdp, (3) 



where Rdp is the ratio of deuterons to protons established during early stages of the fireball. One source of error arises 
due to the neglect of other light composite particles viz. t, ^He and a. Bertsch and Cugnon 2j| proposed to take 
into account these lighter fragments as well generalizing Eq.Q as: 



Sn = 3.945 - In x, (4) 



where 



Hik: 



■JPUke 



d+U+% ■^He + ia 



2" ' 2 



(5) 



p + d + t + 2 3iJe + 2a 

As seen, quantity 'x' measures the yield ratio of deuteron-like (duke) to proton-like (puke) fragments. It has been 
well established from experiments that highest proton multiplicity accounts for most of the charges in HI system. 



thus leaving no room for heavier clusters. Following Ref. 25|, we define the yield ratio of deuteron-like {duke) to 
proton- like (puke) clusters in the following way: 

where Y(A) stands for the number of fragments with mass 'A' in one event. Analogous to experimental results, we 
calculate the total participant multiplicity Np as: 

Np^ f^[Y(A = l) + 2Y{A = 2) 

+3F(A = 3) + 4y(A = 4)], (7) 

where Zp + Zt and Ap + At define the total charge and mass of the colliding system. This procedure allows us 
to calculate the entropy produced in a reaction. We have calculated these ratios within QMD model for unfiltercd 



events, employing minimum spanning tree (MST) procedure of clusterization [26[. In this clusterization approach, we 
assume that two nucleons sharing the same cluster if their centroids are closer than a spatial radius re = |i"a — r/3|. 
One generally choose the radius re = 4: fm. It may be mentioned that some other sophisticated algorithms have also 



been developed which are based on relative momentum of nucleons [2T], fragments' binding energy minimization 28| . 
and backtracing procedure 29] to study the fireball and spectator matter physics. 

III. RESULTS AND DISCUSSION 

Since it is well established that the production of light charged particles and clusters, and ultimately the entropy is 
related to the fireball, it is of interest to see their phase space. In Fig. 1, we display the space (X-Z) of ^^Nb +^^ Nb 
collision at 650 AMeV and at a reduced impact parameter b/b,nax—0-Q- The chosen time represents the freeze out time 
when density reaches asymptotic value. First of all, consistent with earlier attempts, we see that heavier fragments 
belong to the residue of either projectile or target, whereas lighter clusters like free nucleons and light charged particles 
LCPs [2 < ^ < 4] are produced due to the coalescence and emerge from the mid-rapidity region. As shown in Ref. 
[20J , these light charged particles carry vital information about the stopping as well as thermalization of the nuclear 
matter, therefore, are also good candidates for the study of production of entropy in HI reactions. A very little 
influence can be seen of momentum dependent interactions. 

Next we study the final state composite particle yield ratios X/p for the soft (S) and soft momentum dependent 
(SM) interactions. This is shown in Fig. 2 for the collisions of ^^Nb +^^ Nb at incident energy of 650 AMeV as a 
function of impact parameter. Here 'X' stands for A=2, 3 and 4 clusters. On can see that curves for different X/p 
ratios exhibit similar trends: 

(i) The X/p ratio decreases with impact parameter (alternately, increases with Np) indicating more production in 
central collisions compared to peripheral collisions. As shown by many authors [5J, |6|, |25J , Np remains same for nearly 
central collisions and decrease sharply for semi-central collisions and peripheral collisions. 

(ii) For central impact parameters (or, higher Np values), X/p ratios reach an asymptotic value indicating that for 
central events, small variation in impact parameter does not give different results. 

(iii) Role of momentum dependent interactions is nearly marginal justifying the earlier attempts |3(J] and use of soft 
equation of state. 

Since entropy production is mostly measured for central collisions, the use of momentum dependent interactions 
will not give different results compared to static soft equation of state. These different yield ratios implied that one 
also obtains different behavior of density reached, collision rate and multiplicity of various light mass fragments for S 
and SM interactions. 



In Fig. 3, we extend the above study by including the ratio of deuteron-hke (dnfee) to proton-hke {/puke) clusters. 
The calculations for duke/puke are made using soft (S) equation of state for the collisions of ^°Ca +^° Ca (at 400 and 
1050 AMeV) and ^^Nb +^^ Nb (at 400 and 650 AMeV) as a function of participant proton multiplicity. The results 
of the Plastic Ball experiments [6| are also displayed for comparison. The Plastic ball data takes into account the 
overlap region for the yield of deuteronlike and protonlike clusters, while our ratios are calculated for the unfiltered 
events using MST procedure. The yield ratios are calculated typically after 40 fm/c, when average nucleonic density 
saturates and n-n collisions practically cease. At this time, yields of composite particles is well established and may be 
compared with experimental data. One clearly sees that our model describes well the functional form of experimental 
diike/piike ratio which is found to increase with Np (or centrality of the collision) and saturates at higher multiplicity 
end. At low Np, there is a large drop in the yield ratio for model calculations as also found for individual cluster- 
to-proton (X/p) ratios (See Fig. 2). These trends are closely related with nuclear matter stopping and flow effect in 
the formation of hot and dense fireball. Recently, Dhawan et al. [20| studied impact parameter dependence of light 
charged particles (LCPs) yield and anisotropy ratio. It was found that LCPs production was maximum at central 
collisions where maximum stopping of nuclear matter is also achieved. Thus, production of light clusters can act an 
indicator of global stopping achieved in the nuclear matter. Interestingly, duke/puke ratio calculated using dynamical 
approach is in good agreement with experimental data. This shows that one can reliably explore the applicability of 
dynamical approach such as QMD model to further investigate the formation of fireball at intermediate energies. 

It should be emphasized that model of Kapusta [3|, Ij] is highly idealized one that relies on the assumption of thermal 
and chemical equilibrium in the system, for extraction of Sn using Eqs. ([3]) and (|3]). The heavy-ion collisions at higher 
incident energies evolve from an initial non-equilibrium momentum distribution to a thermalized distribution. In a 



heavy-ion system, full equilibrium is not, however, always possible even at central geometry 31|. At lar ge i mpact 
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It is 



parameters, the momentum distribution of nucleons further deviates from the degree of full equilibrium 
worthwhile to study the equilibration of the nuclear matter produced at the time of fireball formation and role of 
secondary de-excitation of heavier clusters on proton-like and deuteron-like abundances during the expansion stage. 
The anisotropy ratio, which is also a measure of the degree of equilibration reached in a heavy-ion reaction, is defined 



as 



201 



{Ra) = F^ . (8) 

The full equilibrium corresponds to {Ra) value close to unity. In Fig. 4(a), we plot the time evolution of anisotropy 
ratio Ra for the reactions of Nb+Nb at incident energies of 400 and 650 AMeV, respectively. The shaded area 
corresponds to the time zone for highly excited nuclear matter which is followed by the decompression phase. Beyond 
this region, the nucleon density saturates and hard n-n collisions cease almost. One can see that {Ra) values approach 
close to asymptotic value just after the compression (t^ 40-45 fm/c) irrespective of the incident energy chosen. The 
{Ra) values at both incident energies are still less than 1. It clearly shows that the HI system has still not achieved 
full equilibrium at the expansion stage of the collision. It also indicates that n-n collisions occurring at later times 
don't alter the momentum space of nucleons significantly. In the lower panels (b)-(c), we display the yields oipuke and 
diike clusters. One can see that in low-density phase, the secondary emission from heavier clusters tend to populate 
the abundances of deuteron-like and proton-like clusters slightly. 

Since full thermodynamical equilibrium is not visible for the participant zone, we shall estimate apparent entropy 
^Sapp^ produced in the HI reactions using Eq.Q. We display in Fig. 5, the apparent entropy Sapp as a function 
of participant proton multiplicity Np for the reactions of Ca 4- Ca, Nb + Nb, and Au -I- Au at incident energies of 
400 and 650 AMeV. One can clearly see that at a given beam energy, it is the volume of participant nucleons (that 
is, Np) which governs the apparent entropy produced rather than the total number of nucleons in the phase space. 
This clearly brings out the participant-spectator picture of HI collisions at relativistic beam energies. It depicts that 



participant volume is solely contributing towards deuteronlike and protonlike cluster abundances independent of the 
target-projectile combination. These results are in agreement with the experimental data and theoretical approaches. 
It may be mentioned that apparent entropy of the fireball would approach realistic baryonic entropy using Eq.(j31) only 
when full thermodynamical equilibrium is achieved for the participant zone. 

Finally we display in Fig. 6, the apparent entropy 'Sapp^ extracted for the central collisions of Ca + Ca (at 400 
and 1050 AMeV) and Nb + Nb (at 400 and 650 AMeV). Also shown here is the baryonic entropy obtained in Plastic 
Ball experiment |6| using the model of Kapusta. Our model entropy values are in close agreement with experimental 
observations ISl, |32l- Nearly no effect of beam energy is visible in experimental as well as in model entropies. 

IV. SUMMARY 

In this paper, we have modeled the dynamical emission of light clusters in highly excited nuclear matter and 
estimated specific entropy based upon the formalism proposed by Kapusta et al. Our model calculations are done 
within the framework of quantum molecular dynam,ics (QMD) model for unfiltered events. The comparison of model 
calculations of yield ratios as a function of participant proton multiplicity Np with experimental data gives quite 
encouraging results. From these results, we conclude that it is the volume of participants in a reaction system that 
governs deuteronlike and protonlike cluster abundances in early interaction stage. The time evolution of momentum 
anisotropy ratio shows that full thermodynamical equilibrium isn't, however, achieved even at later expansion stage 
of the heavy-ion collisions. Our model calculations for apparent entropy produced in the central collisions of Ca+Ca 
and Nb+Nb indicate close agreement with experimental results based upon Kapusta's formalism. Nearly no effect of 
beam energy is visible in experimental as well as in model entropies. 
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Figure Captions 

FIG. 1. The snapshots of spatial coordinates of nucleons in X-Z plane for the reaction of Nb(650 AMeV)+Nb at 
reduced impact parameter b/bmax= 0.6. Time taken here corresponds to the case when n-n collisions cease almost 
just after the violent phase of the reaction. Different colors depict the free particles, nucleons bound in light charged 
particles as well as in heavier fragments. 

FIG. 2. The yield ratio of light clusters to protons (X/p) as a function of impact parameter b using minimum 
spanning tree procedure. The term 'global' for ratios X/p signifies that particle yield is calculated taking full 
ensemble into account and not the limited region only. 

FIG. 3. The duke/puke ratio as a function of baryon charge multiplicity Np. The model calculations (open 
symbols) at the time of freeze out are compared with experimental data (solid symbols). The results are shown here 
for the reactions of Ca+Ca (l.h.s.) and Nb+Nb (r.h.s.). 

FIG. 4. The anisotropy ratio (Ra), and the yields of deuteron-like {duke) and proton-like {puke) clusters obtained 
during the evolution of central ^^Nb +^^ Nb collisions at incident energies of 400 and 650 AMeV. 

FIG. 5. Apparent entropy per nucleon Sapp as a function of baryon charge multiplicity Np for the reactions of 
'^^Ca -1-'*° Ca (open circles), ^^Nb +^^ Nb (half filled circles), and ^^''^Au +1^'^ Au (open squares). Calculations shown 
here are at incident energies of 400 (top) and 650 (bottom) AMeV. 

FIG. 6. The relationship between baryonic entropy and beam energy Eiab for the central collisions of ^"Ca-I-^'^ Ca 
and ^^Nb +^^ Nb. Also shown are entropy values extracted by the Plastic Ball group 6] using Kapusta's formalism. 
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